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Single-particle removal reactions are becoming an important tool for 

studying radioactive nuclei. The nuclei far from stability may reveal new 
regions of large deformation. We discuss the influence of the projectile de- 
formation on the stripping cross sections and longitudinal momentum distri- 
butions of the core. The Glauber formalism for describing such reactions is 
developed and compared with the simpler geometric approach and with ex- 
isting experimental data. The most significant effects are expected for nuclei 
with valence spherical orbitals strongly mixed by deformation. The shape of 
the momentum distribution, in turn, can be used for determining the degree 
of deformation of the projectile. 



I. INTRODUCTION 

Particle removal reactions of various types (nucleon transfer, knockout, stripping and so 
on) have for a long time been one of the most powerful tools for investigation of nuclear 
structure Jl],0]. We are here concerned with a class of reactions at high energies ("stripping" 
in the sense of the word used originally by Serber ||) of which the deuteron breakup is the 
archetype. The strong current interest in these reactions is that they, in inverse kinematics, 
provide a unique opportunity to obtain simultaneously partial probabilities for populating 
individual final states leading to spectroscopic factors as well as momentum distributions 
characteristic of the quantum numbers of the removed nucleon. This allows one to access 
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detailed information on the internal structure of projectile and residue. 

Nuclear breakup reactions are widely used nowadays for studying radioactive nuclei 
@-§H| . A new method was suggested recently for the spectroscopy with radioactive beams 
[/J which is well suited both for "normal" and for halo nuclei. The method uses the longitu- 
dinal momentum distribution of the projectile residue ("core") for the determination of the 
orbital angular momentum of the stripped nucleon. The identification is possible via the 
comparison of the experimental distribution with the theoretical one calculated, for instance, 
in the framework of the Glauber model f9UH|. The association of the definite orbital mo- 



mentum value with the shape of the longitudinal momentum distribution is clearly justified 
for spherical nuclei. First, for a spherically symmetric unpolarized projectile, the role of the 
spin-orbit interaction is not essential for the stripping mechanism. Second, the longitudinal 
momentum distribution is sensitive to the projection m of the orbital momentum of the 
stripped particle onto the beam direction. As a rule, the largest of the possible m-values 
for a given I value dominates, and this I- and m-dependence allows one to determine the 
orbital momentum more or less uniquely |TT |. The results of a recent series of experiments 



with spherical or weakly deformed nuclei |7| confirm the applicability and flexibility of the 
method. 

Recently the first attempt has been made to apply the same techniques to strongly 
deformed projectiles. The inclusive reaction 9 Be( 25 Al, 24 Mg)X was investigated in 0. Many 
neutron-rich nuclei, for example neon and magnesium isotopes in the region N « 30 (the 
island of inversion fl2|), are expected to be strongly deformed, and the role of intruder 
configurations, such as the [330 1/2] level coming from the p/-shell, is extremely important 
13fl . This follows from the theoretical calculations (|14| , |T5| and from the results of the latest 



measurements of the transition probabilities B(E2; 0+ s — > 2 + ) in neutron-rich 26,28 Ne and 
30,32 ' 34 Mg isotopes via intermediate-energy Coulomb excitation ||16|| . Available information 
on the structure of the nuclei in the inversion island is scanty, and one can hope to extract the 
shape parameters from the analysis of the stripping reactions and corresponding momentum 
distributions. 



However, in order to correctly interpret experimental data for the deformed nuclei, one 
needs to modify essentially existing theoretical approaches. Several features of nuclear de- 
formation should be taken into account consistently in the reaction description: (i) the 
deformed mean field of the core is responsible for the strong mixing of different orbital 
momenta of the stripped valence nucleon; (ii) although in the case of axial symmetry the 
projection of the nucleon total angular momentum on the symmetry axis is still a good 
quantum number (usually denoted as Q), the energy splitting and therefore difference in 
the population of different doubly degenerate (±f2) levels preclude complete averaging and 
stress the effects of the spin-orbit coupling; (iii) pure intrinsic f2-states of the nucleon have 
no definite angular momentum projection in the laboratory reference frame defined by the 
beam direction. All mentioned elements are well known but we have to incorporate them 
consistently in the reaction formalism. The realistic calculations along these lines imply 
cumbersome and time consuming numerical work. This might be one of the reasons why, 
in spite of the growing interest in reactions with deformed nuclei [|lj],[n|], our attempt is 
apparently the first one of this kind. 

Of course, scattering by nonspherical nuclei has been studied long ago. The most general 
formalism for this purpose is the coupled channel method outlined in the book [[H]]. Unfor- 
tunately, it is rather hard to use this technically complicated method for the description of 
realistic reactions. A considerably simpler approach, namely the adiabatic approximation, 
was proposed by Drozdov and Inopin in the middle of the fifties pO| , pT| . They used the par- 



tial wave expansion of the scattering wave function and assumed that the collective motion 
(rotation or vibration) of the deformed target nucleus is frozen during the collision time, 
being very slow compared to other degrees of freedom. The term "adiabatic" refers to the 
orientational and shape degrees of freedom only and should not be confused with the adi- 
abaticity with respect to valence nucleon motion. The latter assumption is also frequently 
used in Glauber theory for the description of the reaction mechanism. The Glauber adia- 
baticity means that, starting at energies of 50-100 MeV per nucleon, the characteristic time 
scales for the relative projectile-target motion and for the nucleon internal motion inside the 



projectile are getting very different. 

In the framework of the orientational adiabatic approximation, it is impossible to exam- 
ine the deep excitations of the core. At the same time, soft excitations, with no essential 
reorganization of the internal core structure, are quite treatable. Typical energies of col- 
lective rotational states are low |TS[ compared to characteristic reaction energies. In some 
sense we have here nearly the same situation as in perturbation theory for close or degenerate 
states. To be consistent with the adiabatic approximation, we fix the orientation of inter- 
acting objects during the collision. The state with a certain orientation is a wave packet of 
stationary angular momentum states. This means that after the collision various rotational 
members of the superposition can be excited. 



At the end of the seventies, Abgrall et al. [22] combined the coupled channel formalism 



with the Glauber scattering theory to describe the excitation of the ground state rotational 
band by medium energy (1 GeV) protons. The regular treatment of the adiabatic approxi- 
mation in the framework of the Glauber theory was given in 1990 by Faldt and Glauber ||23| . 



They treated only soft rotational excitations without restructuring of the internal projectile 
wave function. The deformation effect on the stripping cross section was studied only for 
the relativistic deuteron stripping (see also for a qualitative discussion of the defor- 



mation effect on the tensor analyzing power for light projectiles). The approach of our work, 
being a natural continuation of the guideline [E3], is, on the other hand, close to the for- 



mulation suggested in the recent paper by Christley and Tostevin |18| . These authors came 
to a conclusion that the deformation is not very important for the total (angle-integrated) 
reaction cross sections: the corresponding change is less than 1%. Our results are consis- 
tent with this conclusion. However, the longitudinal momentum distributions and partial 
cross sections for stripping from deformed orbitals, which are a strong mixture of spherical 
single-particle states, are sensitive to deformation. 

The paper is organized as follows. The semi-microscopic reaction formalism, taking into 
account the deformation of the projectile, is developed in Sec. ||; the influence of the de- 
formation on the stripping reaction is analyzed, and expressions for the cross section and 



the longitudinal momentum distribution of the projectile residue are obtained. Sec. |Tl| is 
devoted to the study of the deformation effects on the core longitudinal momentum distri- 
bution in the inclusive, with respect to the rotational states of the core, formulation. This 
study is performed using a geometric approach || which is simpler and more transparent 
than the rigorous treatment of Sec. |TI[ In Sec. [TV] numerical calculations of the stripping 
cross sections and longitudinal momentum distributions are presented and discussed in the 
framework of the approaches of Sec. |T| and Sec. HJ for different final states of the residual 
nucleus in a typical reaction 9 Be( 25 Al, 24 Mg)X. The interplay of the internal structure of the 
projectile and the reaction mechanism is illustrated for the model 9 Be( 28 P, 27 Si)X reaction. 
We discuss the perspectives for further studies and propose new experiments in a possible 
region of large deformation in neutron-rich Mg isotopes. Summary and final remarks are 
made in the concluding section. 

II. STRIPPING CROSS SECTION FOR A DEFORMED PROJECTILE 

In the description of projectile breakup reactions, the established terminology is slightly 
different from that conventionally used in the theory of direct reactions Two processes 
are usually distinguished, stripping (or absorption) and diffraction (or elastic breakup). 
In this word usage (see, for example, p6||27| ), elastic breakup is the process in which the 
projectile core and the valence, in particular halo, particle are in the continuum states 
after the interaction with the target. In addition, the projectile residue (core) and the 
target remain in their respective ground states. We generalize this definition to include 
possible "soft" excitations of the core nucleus, such as collective rotations, which disturb 
the intrinsic structure in a minimal way. The absorption differs from the elastic breakup by 
possible excitations of the core and disappearance of the valence particle in the final state 
due to the absorption by the target. In this work we restrict overselves by the treatment 
of the stripping/absorption type of breakup reactions only. This is, generally, the dominant 
channel ||], except for the halo states. Before considering the reaction in detail, let us agree 
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on a general framework of the approach. 

Target and core nuclei are treated as structureless particles in the actual calculations. All 
the observables depend on the internal structure of the core and the target only indirectly 
via corresponding interaction potentials and density functions. Nevertheless, we retain the 
explicit dependence on the internal Jacobi coordinates of the colliding nuclei in the derivation 
of the basic formulas. This makes the derivation more transparent and shows the way for 
extending the approach. In addition, the core may be deformed, in contrast to previous 
treatments p6| , P8|,[l0i . Therefore we need to fix its orientation (more precisely, the orientation 
of the core mean field which determines the geometry of the orbitals for the valence nucleon). 
Although we assume here that the core is axially symmetric, the approach can be generalized 
for triaxial deformation that would be more suitable for nuclei such as, for example, 24 Mg. 
In a similar way, the deformation of the target can be accounted for as would be necessary 
for transfer reactions. 

The interaction of a valence nucleon (for definiteness, a proton) with the target, p+T, 
and the core-target, C+T, interaction are described in terms of phenomenological optical 
potentials. For the C+T interaction, the convolution of a nucleon-target phenomenological 
potential with the core density is used. We do not employ the traditional Glauber approach 
based on the nucleon-nucleon scattering amplitudes |29| , |30| because of relatively small kinetic 
energy per nucleon in the projectile in the experiments we try to describe. It is known, 
however, that the region of validity of the Glauber approximation extends considerably 
beyond what one could expect from the approximations made in the standard derivation 
|T0| . Therefore, our approach can turn out to be very similar to that of Glauber. 

The calculations are performed either in the "laboratory" frame, with the origin coin- 
ciding with the projectile center-of-mass O and z-axis along the beam direction, or in the 
body-fixed frame, characterized by the origin at the core O' and z'-axis along the core sym- 
metry axis. The positions of the target in these frames are given by the vectors R and — r c , 
respectively. The set of coordinates used in our calculations is shown in Fig. [TJ. 
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FIG. 1. Coordinates used for the description of the reaction. 



The description of stripping reactions with exotic nuclei has its roots in the early works 
on the deuteron stripping PJTLf. Attempts to construct a more regular theory of stripping 
processes have appeared considerably later []27j|32| , |33| . Referring to the deuteron stripping 
reaction formalism as a suitable framework for exotic nuclei fragmentation |33| , it is assumed 
by analogy that a stripped projectile constituent (proton) interacts strongly with the target 
nucleus. The core of the projectile escapes as an analog of the neutron. The only difference 
is the possibility for the core to be excited during the collision. The total stripping cross 
section in such an approach is analogous to the proton stripping from the deuteron ||. 
Below we sketch briefly the formal derivation of the total stripping cross section and the 
core longitudinal momentum distribution. 

The projectile (core + valence nucleon) initial wave function is 



7 



|Xi>=c- <Ko - R |0i>, (1) 

where the plane wave describes motion of the center-of-mass with respect to the target, and 
the initial intrinsic wave function of the projectile 

|0 i > = E4°W o ' )t |O> (2) 

VOL 

is a sum of antisymmetrized products of the core internal wave functions generated by the 
operator $wt anc l the function of the relative proton-core motion (creation operator a\,). 
Here the superscript a enumerates various core states, \u) are single-particle orbitals, and 
Affi are corresponding spectroscopic amplitudes. The calculation of the spectroscopic am- 
plitudes for all possible combinations {vol) requires very specific and hard work. Frequently 
only one possibility is allowed or is taken into account, which corresponds to a definite con- 
figuration, for example to the ground state of the spherical (deformed) core + the proton 
occupying a certain spherical (deformed) orbital. This approximation works fairly well in 
the case of light or halo projectiles when there are no bound excited states. The situation is 
more complicated for projectile nuclei which possess many particle-bound states. We limit 
ourselves by using the same approximation with the spectroscopic amplitudes calculated in 
the framework of the particle-rotor model, consistent with our general approach, or bor- 
rowed from the shell-model calculations []3q| . In the case of a deformed core, the internal 
wave function contains information about the orientation of the core symmetry axis with 
respect to the space-fixed frame (the corresponding spherical harmonic Y^Ml- {&a)), and the 
single-particle quantum states \v) are quantized onto the intrinsic symmetry axis. 

The situation with the final state wave function is even more complicated because here 
we have to sum over all possible final configurations of the core and valence proton, including 
the ones where the knocked out proton is absorbed by the target. We do not consider the 
exact structure of this wave function. Instead we determine @,[n|] the absorption cross 



section through a deficit of the probability for the proton to be found at some point r p that 
appears as a result of the p+T interaction. Such an approach is based essentially on the 
spectator mechanism of the stripping reaction and is adequate only in a narrow kinematic 
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region. This mechanism presupposes that the proton is stripped off the projectile while the 
remaining core fragment continues its motion essentially without being disturbed except 
for a smooth action of the optical core-target potential Vc+t taken into account by the 
eikonal core phase shift. The incident energy of the projectile has to be sufficiently large 
and the scattering angle of the core must be small to ensure the applicability of the eikonal 



approximation |27| , |32| . It is assumed in the framework of the spectator mechanism that the 
nucleus-nucleus interaction time is too short for a change of the core symmetry axis direction. 
As we mentioned, this means the possibility to reveal various rotational components of the 
frozen orientation. Therefore the final-state deformed core wave function $W\0) contains a 
spherical harmonic YljM l (^a)j which defines the same direction of the core symmetry axis 
as before the projectile breakup (adiabatic approximation), but in general with Lf,Mf ^ 

According to the idea of the Glauber approximation and our hypotheses of the reaction 
mechanism, the reaction yield is determined by the profile functions for the propagation of 
the valence proton, S p , and the core, S c , in the field of the target potentials, p+T and C+T, 
respectively pBfl . 



S p (b p , z p ) = exp ( -- f P V p+T {h p , z')dz'A , (3) 

V V J— oo / 

S c (b c , z c ) = exp ( -- / Vc +T (bc, z' c )dz' c ) . (4) 

\ V J —oo / 

Here and below we have adopted units with % = 1. The coordinates r p and r c are expressed 
from Fig. [I] as 

A c 1 
r P = "i 7 r it. r c = -r - R. 

P A c + 1 A c + 1 

The impact parameter vectors b p and b c are their corresponding two-dimensional projections 
onto the plane perpendicular to the beam direction. 

In the same spirit, the (core + proton) Glauber wave function in the final state, projected 
onto definite states of core and proton motion with momenta k c and k p , respectively, is 



| X ^> = 5 p (b P ,^)5 c (b c ,^)e^- rc e^^$Mt|o) ; (5) 

where 

Sp(b p , z p ) = exp | ^ jf Vp +T (b p , z' p )dz' \ , (6) 



and 



5 c (b c , z c ) = exp Q £ V c +T(b c , z' c )dz'}j . (7) 
Here we again take into account the interactions with the target by means of the eikonal 



phase shifts [32 



The stripping reaction is defined essentially as disappearance of the proton from the 
projectile as a result of the collision with the target. We need to compare the probabilities 
for the proton to be found at some point r p with and without the p+T interaction. The 
probability amplitude for the core in a given configuration a to have the momentum k c and 
for the valence proton, emerging from the intrinsic state \ v) with the spectroscopic amplitude 
A$?\ to be found at the position r p is 

a(k c ,r p ) = A v S p e~' 1 ^ J d\ c e- ik ^ +lK °- R M^(^uS c ^l\0) , (8) 

where r = r p — r c , (j) v {r) = (r|aj,|0) is the function of the relative proton-core motion 
in the projectile nucleus, and we omit the arguments of the profile functions as well as the 
superscript a. The phase integrals in the profile functions S c and S p in (H) are now taken from 
— 00 to 00 and do not depend on the z-components of the core and proton radius- vectors. The 
intrinsic matrix element in (§) and the profile function S c contain the orientation parameters 
of the deformed core which are adiabatically fixed during the collision. Therefore the whole 
amplitude a(k c , r p ) is still an operator in orientation angles. Thus, apart from the stripping 
probability, it will determine the rotational excitations of the remaining core. 

The core final momentum is a sum of (A C /A)~K (the core part of the initial momentum 
of the projectile as a whole) and the internal momentum q corresponding to relative core- 
proton motion inside the projectile. Due to the large mass difference, A c 3> A p , we can put 
A c m A and obtain for the exponent in the integrand (H) 
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-zk c • r c - iK • R = -zq • r c . 

Eq. (g) thus can be rewritten as 

a(k c ,r p ) = A v S p e~^ J (Pr c e^ T °<f> v {T)(0\il> v S c $t\Q), (9) 
Without the p+T interaction, the expression (|^) reduces to 

a (k c ,r p ) = Ae- !kpT " / d 3 r c e^ r ^(r)(0|<£vS c $t |0). (10) 

The probability for the proton to get absorbed in different spatial points r p , provided that 
the core escapes with the momentum k c , is proportional to the difference of the squares 
of the absolute values of the orientational matrix elements of the operators (|ToD and @ 
between the given rotational states of the core, 

\a (k ci r p )\ 2 fi - |a(k c ,r p )|J. = \A U \ 2 (1 - \S p (b p )\ 2 ) 

2 

(11) 

Here we have explicitly written the arguments of the profile functions and preserved only 
those parts of the core wave functions that define the orientation of the core mean field. In 
the integrand we have the product of four factors: two spherical harmonics defining the core 
rotational state before and after the interaction with the target, the core profile function 
(our operator), and the function of the proton-core relative motion <f) u (r). The dependence 
of 0j,(r) upon the core orientation is considered in detail in the next section. 

The approach above ensures the correct limit of spherical symmetry. Indeed, in this 
case all factors in (|TTD, except the spherical functions, do not depend on the orientation 
angles fl a , and the orthonormalization of the spherical harmonics singles out the diagonal 
transition of the core. Another way to consider the core orientation would be to retain 
the spherical harmonics up to the final calculation of the stripping differential/total cross 
section and then average over all possible directions of the core symmetry axis. This would 
however imply ignoring interference contributions between the processes of rotational recoil 
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J d?r c e-^ J dn a Mr)YZ f M Lf (n a )s c (h c )Y LiMLt (n a ) 



of the core occurring at different positions of the core. In addition, in that case we encounter 
a technical problem of summation over all possible final projections Ml. (final channels). 
Since the values of Mr,, are connected with the orientation of the core axis Q a , the simple 
summation can lead to the double counting. 

The differential cross section for the core final momentum d 3 k c = d 3 q is obtained, as in 
the deuteron stripping theory [§], after integrating the difference of probabilities ([Tl| ) on the 
interaction plane, orthogonal to the projectile direction, over the proton impact parameter 
bp, 

^~EIA| 2 /A(i-l^(b P ) 



(12) 



where J2 means, as usual, the averaging over the initial angular momentum projections and 
summation over the final ones. Here we made a change of variable from the set r p and r c to 
the set r p and r = r c — r p . The three-dimensional momentum distribution of the core (|12| 



differs from the results of refs. |26| , |T0|| by the presence of the spherical harmonics with the 
additional integration over Q a . 

The longitudinal momentum distribution of the core is obtained by integrating over the 
transverse components q^, 



x Jdxdy |^e-^ z |^ a 0,(r)F£ /ML/ (fiJ5 c (bjF LiMii (fi a ) 



The total proton stripping cross section follows after the integration over q z 



(13) 



a = Y}M 2 d%(l-\S p {h p )f) d 3 r 



dn a Mr)Yi fAlLf (n a )s c (b c )Y LiMLi (n.; 



• (14) 



The expression obtained in a number of works f^jTIJ differs by the absence of the angular 



integral. In fact, these results are known from the theory of deuteron stripping [@,|3l[], of 
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course again without deformation effects. The more complicated and formal approaches, 



such as, for example, based on Ref. |[27|| , do not bring any real advantage. 

The presence of the angular integral in ( p!3| , p^|) , together with the angular dependence 
of the core profile function S c (b c ) and, especially, the proton wave function 4> u (r), makes 
a practical calculation of the longitudinal momentum distribution ( |T3"D and the total strip- 
ping cross section ( |14|) a difficult computational problem. Eq. (|i~4f) , for example, contains 



seven embedded integrals. In addition, the function 0„(r) and its arguments depend on the 
angular variables in a very complicated way. To ensure the resonable accuracy one needs 
to carefully check the convergence of each numerical integration. In practice we have used 
a time consuming combination of the grid integration with the Monte Carlo methods. We 
guess that these technical problems were and remain up to now the main obstacle preventing 
the consistent account of the deformation effects in the calculation of reaction observables. 

III. "GEOMETRIC" DESCRIPTION: LONGITUDINAL MOMENTUM 
DISTRIBUTION OF THE CORE 



The longitudinal momentum distribution may be calculated directly from (0). Techni- 
cally, this is very time-consuming. Indeed, we need to perform a multidimensional integra- 
tion over single-particle and core coordinates and over the orientational angles Q a . Prior to 
that, we have to calculate the core profile function S c as a convolution of the nucleon-target 
potential with a core density distribution which, in turn, strongly depends on the orienta- 
tion of the core symmetry axis. To avoid these technical difficulties, one can use a simpler 
approach developed in || on the "physical" level of rigor. 

In the approach |9| we assume that the profile functions are step functions ( "black disc" 
approximation). Further, the interaction time is very short (a sudden approximation) so 
that the core after reaction retains the same momentum distribution as in the projectile 
before the collision with a target. Due to conservation of the total momentum, the core 
momentum distribution is determined by that for the stripped valence nucleon inside the 
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projectile. Therefore, one needs to calculate the the three-dimentional Fourier transform 
of the part of the initial valence nucleon wave function which is cut out by the target. 
Geometrically the overlapping region of the valence nucleon wave function with the target 
corresponds |J] to a cylinder with the axis along the beam direction and the radius defined 
H by the effective radius of the target Rt- 

The approach f{| can be extended straightforwardly in order to take into account the 
projectile deformation. It is easy to include the case of deformed targets also but with 
significant computational complications. In what follows we consider nucleon stripping from 
a deformed projectile on a spherical target. 

The core mean-field deformation can be treated in the framework of the standard Nils- 



son model in the projectile body-fixed frame [p!9| , p8| , p9| . Double degenerate single-particle 
orbitals labeled by the asymptotic quantum numbers [Nn z AQ] are obtained by the diagonal- 
ization of the phenomenological nuclear Hamiltonian in the basis of spherically symmetric 
wave functions 

\[Nn z An}} = J2^ A] m}. (15) 

ij 

Note that the spin-orbit coupling is essential for the correct definition of single-particle 
states. For the sd-shell nuclei we have six double degenerate basis states \ljVL) some of 
which (with the same angular momentum projection Q onto the core symmetry axis) are 
coupled due to the static core deformation. The corresponding single-particle orbitals in the 
laboratory frame result from the rotation T(1Z), 

iKr) = T(K)\[Nn z AQ]) = ^a^ 1 £ {6r}>, (16) 

ij n>=-j 

where D J n , n (K) are elements of the rotation matrix ||40|| . The operator 1Z is defined by the 
Euler angles 61,62, 63 but one of these angles is arbitrary in the case of axial symmetry and 
may be set equal to zero. The remaining two are the polar and azimuthal angles between 
the core symmetry axis and the beam direction, 6\ = <p, 62 = 

The three-dimensional nucleon momentum distribution inside the projectile is the 



squared Fourier transform of the wave function (|16|) 
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i>(q) = J d 3 re~ iq - r ^(r). 



(17) 



For a specific reaction, the momentum distribution corresponds to the spatial integration 
restricted by the region G that, in the plane perpendicular to the beam (z) direction, is a 
disk with the center at the point (x = b min , y = 0, z = 0) and the radius Rt- Here b min is the 
lowest possible value of the impact parameter b m i n = R c + Rt, R c = r A 1 J' s is the effective 
core radius (it is supposed that the interactions with the impact parameter less than b m i n 
lead to violent collisions and cannot be seen in the stripping channel). The longitudinal 
momentum distribution is obtained by integrating \tp\ 2 over the transverse q components, 



dW 

dq z 



2 d 2 q ± 1 f 
—— = — / dxdy 
(ZlT) 6 2tt Jg 



dz e iqzZ tp(x,y,z) 



(18) 



Using eq. fllCf ) we obtain 



dW 

dq z 



1 

2^ Jg 



dxdy 



dze iq * z J2 



a 



[Nn z A] 



E M 



n'=-j 



£(/A, 1/2 Q' - A\jn')R nl (r)Y lA (r)xi/2 



n'-A 



(19) 



Taking into account the orthogonality of the spin functions X1/2W-A, after some algebra 
we obtain the final expression for the longitudinal momentum distribution of the stripped 
nucleon 

dW II,, ^ [Nn z A] 

<j=±l/2 Ij W=-j 



dq z 



— J dxdy 



E «!? A] E 4aW(l 0' + a, 1/2 - a\m 



x / dze^ z R m (Kr)Y m+l7 (Kr) 

-00 



The single-nucleon stripping cross section is obtained by integrating over q z 



rr I dxdy / dz 

' G a=±l/2 J -°° 



E 



a 



[Nn z A] 
ljU 



E D^ n (n)(m' + a,l/2 -a\0) 



Q'=-j 



(20) 



xR m (Kr)Y ia+a (m) 



(21) 
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The longitudinal momentum distribution (p0|) and total stripping cross section (|21|) are 
obtained at a fixed direction of the core symmetry axis. We assumed that the core can 
not change its orientation during the interaction time. But we did not impose any other 
limitations on the core state of motion - the core can be cranked by the process of interaction 
obtaining any angular momentum allowed by selection rules. Since we do not specify the 
state of the core exactly we have here an "inclusive" process with respect to the core final 
states. It means, in turn, that the cross section (^Tj) is to be understood as a sum of partial 
cross sections corresponding to the excitation of the particular core rotational states and 
calculated in the framework of the more regular approach of Sec. |J where all interference 
features were fully accounted for. In distinction to that, in the geometric description the 
stripping processes intitated at different positions of the core are added incoherently and, 
for the unpolarized projectile, we have to average the dW/dq z , eq. (p0|), over all possible 
orientational angles <p, It gives an additional double integration and the factor 1/47T in 
(^y,^). We do not have any averaging in the approach of Sec. [TI| where the integration over 
the orientational angles appears as a part of the calculation of the matrix element over the 
core internal wave functions. An attempt to carry out the averaging leads to uncertainties 
in the summation over core final states (double counting mentioned in Sec. [II]). 



IV. RESULTS AND DISCUSSION 



For practical calculations we use the Nilsson Hamiltonian (with constants from |1| 



expressed in the dimensionless stretched coordinates 3G . The coordinate dependence of the 



basis wave functions in (|T^) upon the rotation operator noticeably lengthens the computer 
calculations. The harmonic oscillator wave functions also have an incorrect asymptotic 
behavior, decaying too fast at large distances. This shortcoming is not very serious for 
tightly bound projectiles (with the binding energies >5 MeV). However for the description 
of weakly bound nuclei, the realistic behavior of single-particle wave functions may be crucial 
so that in this case the results obtained with the aid of the Nilsson Hamiltonian should be 
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considered as preliminary estimates. 

We start with the inclusive reaction 9 Be( 25 Al, 24 Mg)X studied experimentally in [[7|. The 
deformed nature of the 24 Mg nucleus |E| and, correspondingly, 25 A1 projectile, is well es- 
tablished. Our analysis based on the expression (|20| ) can be considerably simplified in the 
case of the observation of the core nucleus in its ground state. Indeed, the valence proton 
in 25 Al occupies the Nilsson orbital [202 5/2]. This orbital is generated by the spherical 
j = 5/2, rrij = 5/2 basis function which is not mixed by the core deformed mean field 
with other spherical orbitals from the srf-shell so that we have only a slight deformation 
dependence through the use of the stretched coordinates. 

The longitudinal momentum distributions of the 24 Mg nucleus in the intrinsic ground 
state (frozen orientation) calculated according to (|20"D for different values of the oscillator 
deformation parameter 5 |19| , |38| , |39l| are shown in Fig. |2]. The theoretical curves are normal- 
ized to the inclusive experimental data, the normalization constants differ for different S. We 
stress again that the experimental data are inclusive and represent the sum of differential 
cross sections to all final states with isospin T=0 and T=l. 
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FIG. 2. Longitudinal momentum distributions of 24 Mg in the stripping reaction 
9 Be( 25 Al, 24 Mg)X with the 24 Mg nucleus in the intr insic ground state. The curves correspond 
to different values of the deformation parameter: 5=0 (solid), 5=0.1 (dotted), 5=0.2 (dashed), 
5=0.3 (long dashed), and 5=0.4 (dash-dotted). The experimental points for the inclusive process, 
together with the error bars, from [7] are shown. The initial beam energy is 65 MeV per nucleon. 

Note that the case 5=0 would not directly correspond to the spherical symmetry of 
the projectile. Indeed, to achieve such correspondence we have to take into account the 
degeneracy of three d 5 / 2 sublevels with Q — - > rrij = ±5/2, ±3/2, and ±1/2. In our case we 
move along the only Nilsson [202 5/2] orbital that corresponds to a single value of \mi\ = 2 



11] and, of course, can not account for other orbitals. To be consistent with the Nilsson 



model, we do not draw the momentum distributions for the mixture of different rrij states. 

The total stripping cross sections obtained from (|20|) by integration over q z are shown in 
Table 1. Although the cross section generally decreases with the growth of the deformation, 
in the region 5 ~ 0.3 — 0.5, where the deformation parameter corresponds to the realistic 
value 5 ~ 0.34 |^1|] for the 24 Mg core, the cross section reaches the minimum and starts to 
increase again. This can be considered as a hint toward the explanation of the difference 
between the theoretical cross section calculated in the framework of spherical nuclei and the 
experimental one. Note again that for a comparison with the case of sperical symmetry we 
need to take into account the stripping from two additional single-particle orbitals [211 3/2] 
and [220 1/2]; the corresponding stripping cross sections at 5 = are equal to 18.6 and 25.28 
mb. The standard averaging over the different orbitals gives approximately 21 mb for the 
cross section corresponding to the absence of the deformation which agrees resonably well 
with the value of ~26 mb, calculated for the spherically symmetric projectile and realistic 
valence proton wave function (with correct asymptotic behavior). Our cross sections should 
always be smaller than their counterparts for the functions with correct asymptotics. 

In the calculation of the longitudinal momentum distributions and stripping cross sec- 
tions we have used the following values for the effective target radius and minimum impact 
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parameter: Rt = 2.948 fm (it corresponds to the cross section a = 273 mb for the free 
proton interaction with 9 Be) and b min = 6.094 fm. 



Table 1. The cross sections, in mb, for the stripping of the [202 5/2] proton from 25 Al to 
the 24 Mg ground state rotational band at different deformation parameters 6. 
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19.2 


18.6 


18.2 


18.1 


18.3 



Let us now compare these results with the more elaborate calculations based on the 
regular theoretical approach of Sec. O. We limit ourselves here by the consideration of the 
realistic 24 Mg deformation parameter 5 ~ 0.34 [4"T] under the assumption of axial symmetry. 
In this approach we can calculate not only the total stripping cross section summed over 
all core collective rotational states, as was done above (Table 1), but also the partial cross 
sections for the transitions into different rotational states of the core. If the results of the two 
rather different approaches are consistent, the inclusive, with respect to rotational states, 
cross section calculated according to the prescription of Sec. |T| has to coincide with the 
sum of all partial cross sections obtained in the framework of the more regular theoretical 
approach of Sec. |I[ The cross sections based on eq. (|i~4"|) for the transitions into three lowest 
levels of the + rotational band of 24 Mg are shown in Table 2. 



Table 2. The cross sections, in mb, for the stripping of [202 5/2] proton from 25 Al for the 
three lowest rotational states of the core. The deformation parameter 5=0.34. 
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The total cross section is equal to 22 mb that agrees with the inclusive cross section »s 18.2 
mb calculated in the framework of the simple approach (see Table 1). The agreement would 
be even closer if we would take into account the value of the spectroscopic factor *4 2 in 
We have estimated the spectroscopic factor using the particle- rotor model for the 25 Al 
nucleus and obtained the total cross section which varies from 20.7 to 21.5 mb. 

As for the shape of the core longitudinal momentum distribution, using the full calcu- 
lation of Sec. HI we obtain a nearly perfect agreement with experimental data as shown in 
Fig. |||. The dashed curve corresponds to the transition into the + core ground state, the 
dashed-dotted one - to the transition into the 2 + state, and the solid curve is the sum of the 
transitions to + , 2 + , and 4 + states. Note that all calculations have been performed without 
any adjustable parameters except for the overall normalization. 
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FIG. 3. Longitudinal momentum distributions of 24 Mg in the stripping reaction 
9 Be( 25 Al, 24 Mg + ) 2+,4+)X. The dashed curve corresponds to the transition into the + core state, 
the dashed-dotted one - to the transition into 2 + state, and the solid curve is the sum of the three 
transitions including the weak 4 + (not shown separately). The experimental points together with 
the error bars from [7] are shown. The initial beam energy is 65 MeV per nucleon. 

Fig. ^ and Table 1 seemed to indicate that the deformation is not very important both for 
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stripping cross sections and longitudinal momentum distributions. This agrees also with the 
results of ref. ||18|| . However, as we see from Fig. |3] and Table 2, such a conclusion would be 
premature. The details of the longitudinal momentum distributions can be explained only 
after taking into account the deformation of the projectile. This is clearly seen in Fig. |] 
where we compare the longitudinal momentum distributions calculated in the framework of 
exact theory of Sec. f|, the solid line, in the extended geometric approach of Sec. |T|, the 
dotted line, and in the original geometric approach of Ref. || developed for the spherical 
projectiles, the dashed line. The comparison of three approaches made in Fig. 4 with the 
absolute normalization to the experimental peak is incomplete since the spherical geometric 
approach, as a rule, leads to the cross sections of considerably reduced magnitude (up to 
25%). 

Three features are worth pointing out in connection with the comparison in Fig. f|. The 
first is that although, as noted earlier in our paper, the use of harmonic-oscillator wave 
functions leads to a wrong asymptotic behavior at large distances, this problem should be a 
minor one here. The reason for this is the relatively large separation energy for the proton 
combined with the effect of the Coulomb barrier, both making the tail of nucleon wave 
function unimportant. Therefore the comparison with the calculation |7J, which used a 
wave function from a spherical Woods-Saxon potential with Coulomb interaction included, 
can be expected to be quantitatively meaningful. In this connection we note, secondly, that 
our best approximation appears to account for the symmetric "shoulders" appearing in the 
experimental data. We take this as the first evidence of the necessity of incorporating the 
effects of deformation in the theory of the high-energy stripping cross sections. And, finally, 
as a word of caution, we remind the reader that the experimental data represent a much 
more complex situation: they are inclusive and an average over intrinsic and rotational 
excitations with isospin 1 and 2 and up to 10 MeV excitation energy. 

Below we show that there are other cases when the effects of deformation are more 
significant. 
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FIG. 4. Longitudinal momentum distributions of 24 Mg in the stripping reaction 
9 Be( 25 Al, 24 Mg)X calculated in the framework of three different models: (i) the exact theory of 
Sec. I, solid curve; (ii) the extended geometric approach of Sec. II, dotted curve; (hi) the geometric 
approach for spherical projectiles of Ref. [9], dashed line. The experimental points together with 
the error bars from [7] are shown. The initial beam energy is 65 MeV per nucleon. 

Let us take as a projectile a nucleus, such as 28 P, which has the valence proton in the 
[211 1/2] Nilsson orbital. In reality, the 28 P nucleus is spherical or very weakly deformed 
but we use this example here to illustrate possible deformation effects. The longitudinal 
momentum distributions for 27 Si in the ground state calculated with eq. (p0|) are shown in 
Fig. [5] for five values of the deformation parameter 5. 
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FIG. 5. Longitudinal momentum distributions of 27 Si in the stripping reaction 9 Be( 28 P, Si)X 
with the Si nucleus in the intrinsic ground state. The deformation values are 5 = for the 
oscillator single-particle wave function (wide solid curve), 5=0.1 (dotted), 5=0.2 (dashed), 5=0.3 
(long dashed), and 5=0.4 (dash-dotted). The fat solid curve corresponds to a spherical projectile. 
It is calculated [7] for a Woods-Saxon potential including the Coulomb interaction which gives the 
correct asymptotics of the tail of the nucleon wave function. All theoretical curves are normalized 
to the experimental data which together with error bars are taken from ref. [7]. The experimental 
data in this case have been corrected for contributions to excited states and represent the branch 
to the 27 Si(5/2 + ) ground state alone. The corresponding partial cross section is 21±5 mb [7]. The 
initial beam energy is 65 MeV per nucleon. 

It is clearly seen that the width of the longitudinal momentum distributions increases 
with increasing projectile deformation. The origin of this behavior is the redistribution of 
the strength of different stretched spherical components in the static field of the deformed 
core, see Fig. ^|. 
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FIG. 6. The total and partial longitudinal momentum distributions of 27 Si in the stripping 
reaction 9 Be( 28 P, 27 Si)X with the Si nucleus in the intrinsic ground state, as a function of the 
assumed deformation values: 5 = for the upper panel, 5=0.2 for the middle panel, and <5=0.4 for 
the lower panel. The longitudinal momentum distributions are shown for the proton from different 
spherical orbitals, d 5 / 2 (dotted), Sjyj (dashed), and d 3 / 2 (dash-dotted); solid curves correspond to 
the total distribution. At the absence of the deformation all curves coincide (upper panel). All 
data are shown in the projectile rest frame, so that the momentum scale is different from that in 
Figs. 2-5. 

At small deformations, Fig. ||(b), the si/ 2 -wave is dominant but as the deformation in- 
creases, Fig. H(c), the and especially G? 3 / 2 -wave overcome, and as a result the shape 
of the momentum distributions undergoes essential changes. It makes it difficult (if at all 
possible) to determine the value of the orbital angular momentum of the stripped nucleon 
from the shape of the core longitudinal momentum distribution. At the same time these 
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changes allow one to estimate the value of the deformation parameter 5. We need to bear in 
mind that the oscillator hamiltonian tends to underestimate the spatial tail of single-particle 
wave functions. This leads to a wider momentum distribution compared with asymptoti- 
cally correct wave functions when we use the normalization to the experimental data. This 
is illustrated by Fig. |5] where the fat solid line corresponds to the asymptotically correct 
spherical (5=0) wave finction. The corresponding total stripping cross section is close to 13 
mb which is considerably less than the experimental value a = 21 ± 5 mb. The ratio of two 
cross sections, calculated for the oscillator wave functions and for the functions with correct 
asymptotic behavior, remains nearly constant in all our calculations. The wider solid line is 
related to the oscillator wave function for the spherical case. The corresponding stripping 
cross sections are shown in Table 3. 



Table 3. Cross sections (in mb) for the stripping of the [211 1/2] proton from 28 P at 
different deformation parameters S calculated with harmonic-oscillator wave functions. 
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FIG. 7. Longitudinal momentum distributions of 2 Si for the hypothetical case of the stripping 
reaction 9 Be( 28 P, 27 Si*)X with the proton knocked out from the [202 5/2] orbital; the calculations 
are based on eq. (20). The deformation parameters are 5=0 (solid curve), (5=0.1 (dotted), 5=0.2 
(dashed), 5=0.3 (long dashed), and 5=0.4 (dash-dotted). The experimental points together with 
errors from [7] are shown. The initial beam energy is 65 MeV per nucleon. 
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FIG. 8. The same as Fig. 7, for the stripping from the [211 3/2] orbital. 
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Table 4. The cross sections (in mb) for the stripping of [202 5/2], [211 3/2], and [220 1/2] 
protons from 28 P at different deformation parameters 5. 
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FIG. 9. The same as Fig. 7, for the stripping from the [220 1/2] orbital. 

One of the attractive features of our simple approach is the possibility to treat the 
stripping reactions with the core excitation in the same manner. For example, for the 
particular reaction 9 Be( 28 P, 27 Si*)X with the detection of the excited projectile residue 27 Si*, 
we have three orbitals in the sd-shell, [202 5/2], [211 3/2], and [220 1/2]. Assuming that 
the measured |7| cross section is originated by the knockout from an individual orbital we 
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obtain the longitudinal momentum distributions shown in Figs. [7|-[|. For the comparison 
with experimental data, a — 70 ± 11 mb, we have to sum the partial cross sections because 
the different core final states were not separated in the experiment 0. The corresponding 
values are shown in the last row of Table 4 including the contribution from the ground 
state orbital [211 1/2] given in Table 3. In the treatment with the use of the spherical 
shell model @ it is important to include the relevant spectroscopic factors. In the simple 
Nilsson model used here, we assume that pairing is weak and that the occupation numbers 
simply take the values (0,1,2). We do not take into account the relevant spectroscopic factors 
which depend on the deformation of the core mean field and can not be merely borrowed 
from the standard shell-model calculations. To make our consideration fully quantitative, we 
would need also (i) to correct the asymptotic behavior of harmonic oscillator wave functions, 
(ii) to take into account feeding from high-lying states ||. Nevertheless, even oversimplified 
estimations of the partial and total cross sections, Table 4, show a relatively weak dependence 
on deformation whereas the shape of the longitudinal momentum distributions, Figs. [7|-^], is 
more sensitive to the deformation parameter. 

Based on those results we may try to predict the shape of the longitudinal momentum 
distributions of the projectile residue in reactions with strongly deformed nuclei. As a 
possible candidate for the role of the deformed projectile we can use neutron-rich 26 ' 28 Ne 
or 30 > 32 > 34 Mg isotopes. The measurements and analysis of recent work [FJ| allow one to 



expect a large deformation of these nuclei. Therefore, it would be interesting to study the 
longitudinal momentum distributions of the projectile residue, for example in the reaction 
9 Be( 30 Mg, 29 Mg 9 s )X. The preliminary estimate based on the measurement of the reduced 
transition probability B(E2; 0+ s . — > 2 + ) gives 5 m 0.46 — 0.52 for the deformation parameter 
of the projectile. It corresponds to the knockout of a neutron either from [321 3/2] or from 
[202 5/2] Nilsson orbitals (here we extrapolate the appropriate Nilsson diagram to large 
deformations). The first orbital comes out from the pf -shell (intruder) whereas the second 
one belongs to the s<i-shell. The orbitals intersect just in the region of 5 indicated by the 
experiment [Q . Thus, we obtain an interesting possibility to determine actual deformation 
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of Mg via the comparison with experimental data and clarify the internal structure of 
this nucleus. Similar arguments are true for other isotopes in this region. To illustrate this 



possibility we have calculated, see Fig. |X0j, the longitudinal momentum distributions for the 
neutron knocked out from [202 5/2] Nilsson orbital, solid curve, and from the intruder orbital 
[321 3/2], dashed curve. Two results differ very strongly both in the magnitude of the cross 
section and, which is even more essential, in shape. The estimate of the corresponding total 
cross sections gives ~ 8 mb for the [202 5/2] and ~ 16 mb for the [321 3/2] orbital. 
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FIG. 10. Longitudinal momentum distributions of 29 Mg in the stripping reaction 
9 Be( 30 Mg, 29 Mg)X with the 29 Mg nucleus in the intrinsic ground state. The neutron is knocked 
out from the [202 5/2] orbital, solid curve, and from the intruder [321 3/2], dashed curve. The 
deformation parameter for both curves is 5=0.5. The distributions are calculated in the projectile 
rest frame. 



V. CONCLUSION 

We investigated the influence of the projectile deformation on stripping cross sections 
and longitudinal momentum distributions of the core. For this purpose we developed an 
appropriate formalism treating the deformed projectile mean field as in the standard Nils- 
son model. In practical calculations, we limited ourselves by the fast processes with the 
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core nucleus as a spectator its orientation being adiabatically fixed during the short collision 
time. The study of reactions like 9 Be( 25 Al, 24 Mg 5 . s )X, 9 Be( 28 P, 27 Si g . s .)X, and 9 Be( 28 P, 27 Si*)X 
shows the following. 

(i) The general influence of the deformation on the stripping cross section and the core 
longitudinal momentum distribution is relatively weak for the knockout of nucleons from or- 
bitals which are not, or are only slightly, mixed by the deformed mean field of the projectile. 
Nevertheless, even in such cases it is possible to reproduce the details of the experimental 
longitudinal momentum distributions only taking into account the deformation effects. 

(ii) The deformation of the projectile is seen the most clearly in the knockout of nucleons 
from the strongly mixed orbitals. In the case of the axial symmetry it corresponds to the 
lowest possible fl. 

(iii) The influence of deformation may be less pronounced in the case of halo systems where 
the wave function of the valence nucleon is extended far outside the core. The choice of n Be 
as a projectile |^2| seems to be a good testing ground for verifying this conclusion. 



(iv) The use of harmonic-oscillator wave functions with a wrong asymptotic behavior at 
large distances is quite meaningful for a projectile with high nucleon separation energy. 

(v) The geometric approach to the description of the stripping reactions is quite applicable 
for an estimation of the reaction cross sections but it is oversimplified for reproducing the lon- 
gitudinal momentum distributions where we need to use a more detailed semi-microscopical 
treatment. 

(vi) The developed in Seq. |I| semi-microscopic approach appears to account for the symmet- 
ric "shoulders" of the longitudinal momentum distributions appearing in the experimental 
data. We take this as the first evidence of the necessity of incorporating the effects of de- 
formation in the theory of the high-energy stripping reactions. 

(vii) The analysis of the shape of the core longitudinal momentum distribution allows one 
to estimate the degree of the projectile deformation. This is necessary for determining the 
internal structure of radioactive nuclei. The region of neutron-rich heavy magnesium iso- 
topes would be a promising field for future experiments. 
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(viii) It would be interesting to demonstrate in further calculations the predicted by a sim- 
pler approach difference in the momentum distributions of a strongly deformed core for the 
processes leading to the individual members of a rotational band. 

(ix) In the case of the deformed projectile it is impossible to neglect the spin-orbit interac- 
tion since it defines the single-particle orbitals in the deformed core mean field. 

(x) An analysis of experimental data ignoring the effects of the deformation can lead to 
significant errors in assigning orbital angular momentum involved in the stripping reaction. 
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